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Abstract 

A unitarized nonrelativistic meson model which is successful for the description of the 
heavy and light vector and pseudoscalar mesons yields, in its extension to the scalar mesons 
but for the same model parameters, a complete nonet below 1 GeV. In the unitarization 
scheme, real and virtual meson-meson decay channels are coupled to the quark-antiquark 
confinement channels. The flavor-dependent harmonic-oscillator confining potential itself 
has bound states e(1.3 GeV), S{1.5 GeV), 5(1.3 GeV), k(1.4 GeV), similar to the results 
of other bound-state qq models. However, the full coupled-channel equations show poles 
at e(0.5 GeV), 5(0.99 GeV), 5(0.97 GeV), k(0.73 GeV). Not only can these pole positions 
be calculated in our model, but also cross sections and phase shifts in the meson-scattering 
channels, which are in reasonable agreement with the available data for vrvr, gir and Kit in 
S'-wave scattering. 
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1 Introduction 



The rich structure in meson-meson scattering at intermediate energies has stimulated many 
theoreticians to fit the existing quark models to the experimental results [1-7]. Especially 5"- 
wave meson-meson scattering shows structures which are very intriguing [8-19]. Detailed phase- 
shift analyses reveal two pronounced scalar mesonic resonances below 1 GeV, namely the S'(975) 
resonance in tttt [8-12, 14] and the 5(980) in rjn [15, 16] S-wme scattering. The other relevant 
resonances that appear nowadays in the tables of particle properties [17] are the e(1300) in tttt 
and «:(1350) in Kn [18, 19] 5'-wave scattering. 

It is well known that these particles cause severe problems if one wants to understand them 
as quark+antiquark (qq) states. For instance, confronted with the S'f/(3)flavor quark model, the 
resonance positions do not fit the quadratic or linear Gell-Mann-Okubo mass relations [4] (see, 
however, [3]). A possible resonance in tttt S'-wave scattering at 600 MeV, which was poorly 
recognized in early analyses [8], disappeared from the tables of particle properties in the seventies. 
Nevertheless, some years later this resonance revived within the bag model, due to a solution 
for the scalar- meson problem presented by Jaffe [6] , who pointed out that these resonances stem 
from qqqq states. The large binding energy which is assumed for such configurations makes the 
low masses required by experiment possible. All kinds of quark configurations [6, 20, 21] and 
gluon-gluon bound states [7] might exist, other than the standard qq for mesons and qqq for 
baryons. This probably no one doubts, but there is no experimental evidence that they should 
couple significantly to hadron-hadron scattering [2,5]. 

To select the e(1300) resonance as the isosinglet partner of the 5(975), rather than the e(600), 
is probably the result of bag- model interference with the analysis of the tttt (S'-wave scattering data, 
because in the bag model, and also in other bound-state hadron models, the lowest J^*^ = 0++ 
isospin-zero qq object fits better with a total mass of about 1.3 GeV [20,22]. 

In this paper we will show that we have no difficulties to explain the scalar mesons within 
our unitarized quark model, and to interpret them as qq states with a meson-meson admixture. 
However, neither the model nor the data exclude poles in the scattering matrix which do not 
appear in the tables but nevertheless might be interpreted as resonances. 



2 



2 The Model 



The unitarized quark model is described in many articles. We will therefore confine ourselves 
to only briefiy discuss the main features here and to give a complete list of references, [23-28]. 
In our treatment, meson- meson scattering processes couple to qq quark configurations or mesons 
via the annihilation and creation of a qq pair out of the vacuum. The reverse coupling describes 
the decay process of a meson or the coupling of a meson to its virtual decay channels. In [23] 
and [24], the explicit form of a multichannel Schrodinger description of such a system is given (see 
also Appendix . Several meson-meson scattering channels are, via the QPC mechanism [29], 
coupled to permanently closed qq channels with the same quantum numbers. It is also shown 
in [23] how to account for relativistic effects, and for the effects of one-gluon exchange in the qq 
channels. 

Scattering matrices, phase shifts, cross sections and wave functions can be calculated from 
the Schrodinger equation by an approximative method [25, 26] which leads to an S-matrix that 
is explicitly analytic in the complex energy plane, and unitary. 

Phase shifts and cross sections can be checked to be in good agreement with the data if 
available. In other cases, the pole positions of the scattering matrix can be compared with the 
bound-state and resonance positions found by experiment. Wave functions may be compared 
with those expected from leptonic decays. 

It has been our observation [23] that the properties of the J^*-^ = 1 and J^*^ = ^ mesonic 
resonances are reasonably well described with a few model parameters: the effective quark masses, 
where the effective up and down masses can be taken equal, one universal harmonic-oscillator 
frequency which describes the confining force in the permanently closed qq channels for all possible 
fiavor configurations, and two or three parameters to describe the coupling of the scattering sector 
to the confinement sector. 

In this investigation, we apply the model to S-wave meson-meson scattering. The quark 
and the antiquark in the permanently closed channel (s) move in relative P- waves, whereas the 
mesons in the scattering channels are in relative S- and Z^-waves. For the e and S we use one 
Schrodinger equation with two permanently closed channels, one for the nn pair and one for the 
ss pair. The mixing occurs in our model quite naturally via the coupling to scattering channels 
which contain strange mesons. We will discuss the results furtheron. 

In the first place, we do not alter the effective quark masses, nor the universal harmonic 
oscillator frequency. The only place where we allow some minor changes, if necessary, is in the 
potential which couples the confinement and the decay sectors. For the vector and pseudoscalar 
mesons, the so-called color splitting could be accounted for by a component of this potential. As 
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a result of our calculations for the scalar mesons, we conclude that in their case other possible 
interactions seem to compensate the effects of color splitting. So we decide to set to zero the 
parameter which regulated the color splitting in the case of the vector and pseudoscalar particles, 
and not to take other contributions into account. The only legitimations of this procedure are 
the facts that the results came out reasonable and that it is not very relevant for the point we 
want to make in the present paper. The other two parameters in the coupling potential remain 
unaltered with respect to the corresponding parameters in the case of vector and pseudoscalar 
mesons. 
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3 Results 



Let us first discuss S'-wave tttt scattering. The lowest bound state of our confining potential for 
J^^ = 0++ qq pairs has a mass of about 1.3 GeV, which is at precisely the same place as the 
ground state of other bound-state meson models. If we turn on the overall coupling constant 
of the transition potential, bound states show up as resonances in vrvr scattering. At the model 
value of the overall coupling constant, which was obtained from the analysis of pseudoscalar and 
vector mesons [23], a pole shows up with a real part of about 1.3 GeV, which accidentally equals 
the above-mentioned bound-state mass. Naively we might expect that one would only find a 
resonating structure in vrvr scattering in that energy domain. However, Fig. ([T]) shows that the 
calculated phase shifts have structures at much lower energies, which indicates that low-lying 
resonance poles have been generated. 

deg 
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0:4 0:6 0;8 GeV 

vrvr invariant mass 

Figure 1: Elastic S-wave vrvr phase shifts. The various sets of data are taken from (0, [9]), 
X, o, <, >, for analyses A, B, C, D, and E, respectively, of [10]), (o, [11]), (*, [12]), and (•, [13]). 
The solid line is our model result. 

We can scan the complex energy plane for these poles in the scattering matrix, finding one 
pole at about 450 MeV with a roughly 250 MeV imaginary part, and another pole at the 5(980) 
position. The imaginary part of the first pole is so large that a simple Breit-Wigner parametriza- 
tion is impossible, and large differences between the "mass" of the resonance and the real part 
of the pole position will occur. How these poles are connected to the harmonic-oscillator bound 
states is a very technical story, which is beyond the scope of this paper; suffice it to state that 
such a connection exists. As we have discussed in [27], these poles are special features of S'-wave 
scattering and do not show up in P- and higher-wave scattering, which explains quite naturally 
why they are not found there. 
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Figure ([T]) also shows that the new structures at low energies are in reasonable agreement 
with the experimental situation. A criticism which might come up if one inspects this figure in 
more detail is that the theoretical phase shifts do not fit the data to a high precision, but only 
roughly follow the experimental slope in the data. However, this is not a fair criticism, since 
we are comparing our calculations with the raw data, with unsubtracted background, which are 
presently the only available data. 

In our model we left out of consideration all possible final-state interactions in the scattering 
channels, like meson exchange, Pomeron exchange, quark interchange, etc.. Moreover, the form 
of the transition potential may be too simple, since, for example, we have taken a local transition 
potential, and it could also have a more complicated r dependence due to more sophisticated 
meson-decay form factors. So our calculations better do not follow the data very accurately. It 
remains, however, a pity that no analysis exists for meson-meson scattering which subtracts the 
known effects and leaves us with the consequences of the remaining interactions, a strategy that 
is nowadays popular in analyzing nucleon-nucleon scattering data [30], because then we could 
really see how good the remaining interactions are accounted for in our approach. From our 
present calculations, we must conclude that final-state interactions will probably alter the phase 
shifts a bit in the region around 600 MeV in order to change the slope of the curve towards the 
data. Note that the phase shifts for low energies are almost completely accounted for by the 
coupling to the permanently closed qq channels. 




r]TT invariant mass 



Figure 2: Elastic 5- wave rjir cross section. The solid line is our model calculation. Data are 
taken from [15]. 

In S-wave rjir scattering, the data are limited to cross sections in the energy domain of the 
6. Here we find that straightforward calculations lead to problems with the position of the 6 
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resonance. We suspect that these difficulties are connected to the U{1) problem. Our strategy in 
this case is discussed below, and the result is depicted in Fig. ([2]), where we shifted the calculated 
cross section by 20 MeV in order to get the peak values of the experimental and the theoretical 
curves on top of each other. 




Kit invariant mass 



Figure 3: Elastic I = \ S'-wave Kn phase shifts. Data are from (©, [4] and [18]) and from 
(•, [31])^. The dashed line is our model calculation (not a fit). 

The results for Kn are depicted in Fig. ([3]). We see there that the phase shifts for low energies 
are rather well produced by the model, but at higher energies only a rough description of the 
data is given: the number of resonances at some energies agrees with the data but the detailed 
structure is not reproduced at all. Also here we have difficulties which are presumably related 
to the U{1) problem. 

In previous investigations [23], we took a nonstrange quark content for the f] meson and a 
strange one for the t]' meson, which is called ideal mixing. The data in the case of scalar mesons 
are, however, more sensitive to the quark contents of the rj^s and we found the following: the 
TjK channel in the isodoublet case is much less coupled than what would follow from our general 
approach as described in [28], because the data show that there is not much inelasticity below 
rj'K. The best result is obtained if the rjK channel is completely decoupled and rj'K enhanced 
to compensate. 

Something similar appears to be necessary in the isoscalar case: if we take for rj the ss and 
for 7]' the nn system, we find the best results for the 5* pole, although the complete (e, 5*) system 
is not very sensitive to these changes. The S, however, is very sensitive to our approach, since 
the lowest threshold is rj-R. In this case we have to reduce the rjiT coupling by a factor 1/6 and 
^Added after publication. 
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to enhance the -q'n couphng to compensate. This would be the case if the SU{3) mixing angle 
between rji and rjs was +11°. The plus sign is puzzling, because Tornqvist [2] claimed that he 
needs the standard —11° [17] in order to fit the data. A more rigorous study on this point is 
in preparation, which shows that part of the problem might stem from our choice of transition 
potential. 

4 Conclusions 

If we allow only minor differences in the transition potential and take some action with respect 
to the f] mesons, we find that the scattering matrices for S'-wave meson-meson scattering agree 
reasonably with the data. An inspection of the complex energy plane shows the following poles 
in the various scattering matrices: 5(994 - 20i MeV), 5(968 - 28i MeV), e(470 - 208i MeV), and 
k{727 — 263i MeV). Preliminary results show that the e and k poles are rather sensitive to the 
transition potential. The e is, however, always somewhere around 500 MeV central value, with 
a large width. The k may vary more and even become slightly heavier than 1 GeV [32], with 
a rather large width. Nevertheless, the main conclusion is that these poles occur as normal qq 
configurations with meson-meson admixtures. 

Our calculations show clearly that there is no need to incorporate qqqq channels in our meson 
model. This indicates that there is no phenomenological reason why genuine qqqq configurations 
should couple strongly to meson-meson channels. 
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A Appendix 



The resonances are described by the set of Schrodinger equations 

|-|^ + ^^^^^ + ME)V{r) - k\E)^ Mr) = , (1) 

which consist of n confined (permanently closed) channels and m free (open or closed scattering) 
channels. So L, /i, k"^ are (n + m) x (n + m) diagonal matrices containing the orbital angular mo- 
menta, the reduced masses, and the momenta of the various channels. The latter two quantities 
are determined from 
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I j 1 dk 

E = y/c^ + mi^ + ^Jk'^ + m2^ and /i = " ^'^^ 

Here, mi and m2 stand for the masses of the quarks in the confined channels or the meson 
masses in the scattering channels. The latter are taken to be the experimental values, while the 
quark masses were determined [23], for the J^^ = O"'^ pseudoscalars and J^'" = 1 vectors, to 
be 

m„ = 406 MeV, m, = 508 MeV, m, = 1562 MeV, and nib = 4724 MeV . (3) 

The nonrelativistic limit of ([2]) is used in the confined channels, and for energies lower than 
the threshold also in the free channels. 

The [n + m) X [n + m) potential matrix reads 

V = { . (4) 

The confining potential Vc is a diagonal nxn matrix containing the mass-dependent harmonic 
oscillators 

[Vc]u = l^^^u;'r' . (5) 

The mxm matrix Vj describes possible final-state interactions. The nxm transition potential 
Knt is taken to be 

where g is the universal coupling constant, tq = p^/ ^J^iUJ the transition radius, and E/D a 
phenomenological factor, with D = mi + m2. The above-mentioned parameters were also fixed 
[23] by the J^*" = 0"^ pseudoscalars and J^*" = 1 vectors to 
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u = 190 MeV, g = 10.4, and po = 0.56 . (7) 

In comparison to the 0-+ and J^'^ = 1"" case, we left out the color term, which was necessary 
for the pseudoscalar/vector splitting. This is, of course, not a very serious intervention, since for 
P-wave qq systems the color splitting is anyhow much smaller than for 5- wave systems [32]. 

The numbers Cij are the relative decay couplings. If we pursue the concepts behind the 
confining potential further to the phase of the pair creation, we are led to the assumption that 
all the interquark forces are harmonic-oscillator forces during this process. This idea provides us 
a scheme in which the possible decay channels and their relative strengths can be calculated [28] . 
Their values are listed in Table [TJ 

The rjn and rjg denote the pure nn and ss states, respectively. Of course, these are not the 
physical rj and t]'. See the Conclusions for a discussion about this point. The tttt (S-wave) 
scattering is described in one system of equations, with two confined channels containing an nn 
and an ss pair. These systems are linked to each other via decay in which strange mesons occur. 
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Decay Products 


Initial Mesons 


Channel 


Spin 


e and S 


5 


K 






nn 


ss 






TTTT 





3/40 








VnVn 





1/40 








VsVs 







1/16 






KK 





1/40 


1/16 


1/24 




PP 





1/40 








PP 


2 


1 /o 

1/2 








UJUJ 





1/120 








ijjijj 


2 


1/6 








K*K* 





1/120 


1/48 


1/72 




K*K* 


2 


1 In 

1/6 


5/12 


5/18 




00 







1/48 








2 




5/12 















1/12 

















pu 

r 









1/36 




pu 


2 






5/9 




ttK 











1/16 













1/48 


n K 


n 








1/24 


pK* 











1/48 


pK* 


2 








5/12 


uK* 











1/144 


ujK* 


2 








5/36 


(f)K* 











1/72 


(t)K* 


2 








5/18 



Table 1: Relative quadratic coupling constants (see text) for the decay process of a scalar 
meson into a pair of pseudoscalar/ vector mesons. 
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